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ABSTRACT 

We present abundance measurements from 26 emission-line- selected galaxies at z ~ 0.6-0. 7. By reaching 
stellar masses as low as 10 8 M©, these observations provide the first measurement of the intermediate-redshift 
mass-metallicity (MZ) relation below 10 9 M©. For the portion of our sample above M > 10 9 M q (8/26 galaxies), 
we find good agreement with previous measurements of the intermediate-redshift MZ relation. Compared to 
the local relation, we measure an evolution that corresponds to a 0.12 dex decrease in oxygen abundances at 
intermediate redshifts. This result confirms the trend that metallicity evolution becomes more significant toward 
lower stellar masses, in keeping with a downsizing scenario where low-mass galaxies evolve onto the local MZ 
relation at later cosmic times. We show that these galaxies follow the local fundamental metallicity relation, where 
objects with higher specific (mass-normalized) star formation rates (SFRs) have lower metallicities. Furthermore, 
we show that the galaxies in our sample lie on an extrapolation of the SFR-M* relation (the star-forming main 
sequence). Leveraging the MZ relation and star-forming main sequence (and combining our data with higher-mass 
measurements from the literature), we test models that assume an equilibrium between mass inflow, outflow, and 
star formation. We find that outflows are required to describe the data. By comparing different outflow prescriptions, 
we show that momentum, driven winds can describe the MZ relation; however, this model underpredicts the amount 
of star formation in low-mass galaxies. This disagreement may indicate that preventive feedback from gas heating 
has been overestimated, or it may signify a more fundamental deviation from the equilibrium assumption. 
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1. INTRODUCTION 

The balance between gaseous inflows, outflows, and star 
formation is a critical frontier in our understanding of galaxy 
evolution. Feedback caused by stellar winds, supernovae, and 
supermassive black holes is often used to explain a variety of 
observations, from luminosity and stellar mass functions to 
the enrichment and reionization of the intergalactic medium 
(IGM). However, a complete physical picture of these feedback 
processes (i.e., Murray et al. 2005; Hopkins et al. 2012) is still 
debated. Continued efforts to provide new observational tests 
are essential. 

The correlation between galaxy stellar masses and gas- 
phase metallicities (the mass-metallicity (MZ) relation) is 
one important probe of star formation feedback. The galactic 
outflows that slow star formation and the inflows that promote 
it can also alter metallicities. On one hand, metal-poor material, 
when accreted onto a galaxy, can lower the metallicity. On the 
other hand, supernova-driven winds may remove metal-enriched 
material from galaxies. It has been known for many years that 
models that ignore inflows and outflows (i.e., closed-boxes) 
fail to reproduce observed abundance patterns in galaxies (van 
den Bergh 1962). Indeed, in recent years, both analytical and 
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numerical models have shown that outflows (and sometimes 
gas accretion) are needed to explain the MZ relation (Tremonti 
et al. 2004; Dalcanton 2007; Brooks et al. 2007; Finlator & 
Dave 2008; Erb 2008; Dave et al. 2011a, 2012; Peeples & 
Shankar 2011; Dayal et al. 2013; Lilly et al. 2013). However, 
to date, studies have not converged on the properties (rates, 
kinematics, metal enrichment, halo mass dependence, and 
redshift evolution) of these gaseous flows. 

One avenue to better understand the physics that governs the 
MZ relation is to study low-mass galaxies (e.g., Lee et al. 2006; 
Zhao et al. 2010; Berg et al. 2012; Wuyts et al. 2012). In these 
systems, galactic winds are especially effective at escaping the 
gravitational potential of their hosts, regulating star formation 
and enabling enrichment of the IGM (Oppenheimer et al. 2009; 
Kirby et al. 2011). Hence, observational constraints on low- 
mass galaxies offer some of the most stringent tests of galaxy 
formation models. 

Outside the local universe, the MZ relation is poorly con- 
strained at stellar masses below 10 9 Mq (and 10 10 Mq for 
Z > 1). While large spectroscopic surveys (Le Fevre et al. 2004; 
Lilly et al. 2007; Newman et al. 2012) have enabled abundance 
measurements of statistical samples out to z ~ 1, these surveys 
are typically limited to R < 24. Hence, these intermediate- 
redshift MZ relations have been derived for M > 10 9 Mq (Lilly 
et al. 2003; Savaglio et al. 2005; Lamareille et al. 2009; Zahid 
et al. 2011; Cresci et al. 2012; Moustakas et al. 2011). Never- 
theless, in the case of the Cosmic Origins Survey (COSMOS), 
the extensive broad and intermediate-band photometry allows 
reliable mass constraints an order of magnitude lower at inter- 
mediate redshifts. Therefore, spectroscopic follow-up of fainter 
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galaxies can significantly extend the intermediate-redshift MZ 
relation. 

In this paper we use an emission-line-selected sample to place 
new constraints on the low-mass end of the MZ relation at 
z ~ 0.6-0. 7. By drawing our sample primarily from the ultra- 
faint emission line objects that we have previously identified 
with blind spectroscopy in the COSMOS field (Martin et al. 
2008; Dressier et al. 201 lb; Henry et al. 2012), we obtain oxygen 
abundances for galaxies with stellar masses of 10 8 M Q < M < 
10 1() M q . In this manner, we provide the first constraints on the 
low-mass evolution of the MZ relation, reaching stellar masses 
that are comparable to the limiting mass of local Sloan Digital 
Sky Survey (SDSS) samples. 

This paper is organized as follows: in Section 2 we describe 
our spectroscopic observations and the COSMOS imaging 
data that we use. In Section 3 we describe our emission-line 
measurements and stellar mass derivations. Then, in Section 4 
we calculate the oxygen abundances of galaxies in our sample, 
discussing the various diagnostics that have been proposed 
to break the degeneracy of the double-valued R23 metallicity 
indicator (Pagel et al. 1979). In Sections 5 and 6 we compare 
our MZ relation to previous derivations and investigate the 
presence of a mass-metallicity-SFR fundamental plane. Finally, 
we compare to theoretical predictions of the MZ relation in 
Section 7. In this paper we use AB magnitudes, a Chabrier 
(2003) initial mass function (IMF), and a ACDM cosmology 
with Q. m = 0.3, £2a = 0.7, and Hq = 70 km s -1 Mpc -1 . 
Throughout the text we report measurements of doublet lines: 
[On] 113727, 3729 and [Oih] 114959, 5007. For the sake of 
brevity, we use the notation “[O n]” and “[O hi]” to refer to both 
lines in the doublet, or, when appropriate, the sum of their fluxes. 

2. OBSERVATIONS 

2.1. Target Selection and Follow-up Spectroscopy 

The emission-line galaxies in the present sample were initially 
identified as part of our multislit narrowband spectroscopic 
survey. The observations are presented in detail in Martin et al. 
(2008), Dressier et al. (2011b), and Henry et al. (2012). In 
brief, our design uses the Inamori-Magellan Areal Camera and 
Spectrograph (IMACS; Dressier et al. 2011a) on the 6.5 m 
Magellan Baade Telescope at Las Campanas Observatory. We 
used a Venetian blind slit mask and narrowband filter centered in 
the 8200 A OH airglow free window. This method allows for the 
efficient selection of emission-line galaxies and reaches fluxes 
as low as 2.5 x 10” 18 erg s“ 1 cm -2 — a factor of five fainter than 
narrowband imaging surveys (e.g., Kashikawa et al. 2011). 

Because of the faintness of the emission lines and the narrow 
bandpass in our search data, we use follow-up spectroscopy 
with DEIMOS (Faber et al. 2003) to identify the redshifts 
of the IMACS-detected galaxies. The primary goal of these 
observations was to confirm Lya-emitting galaxies and measure 
the faint-end slope of the Lya luminosity function at z = 5.7 
(Henry et al. 2012). However, for the foreground galaxies at 
z ~ 0.6-0. 7, these follow-up observations contain the [Om] 
114959, 5007, [Oil] 113727, 3729, and H/3 emission lines that 
comprise the R23 metallicity indicator (Pagel et al. 1979). Three 
classes of these objects were included on follow-up slit-masks. 

1 . Emission-line galaxies for which continuum in the search 
data ruled out the Lya identification were chosen as MZ 
targets. To obtain the best possible spectrum, we searched 
the COSMOS photometric redshift catalog (Ilbert et al. 
2009) for nearby galaxies (<2") with z ~ 0.6-0. 7 and 


shifted the slit positions to coincide with these matches.' 1 In 
total, the final sample of 26 star-forming galaxies (described 
in Section 3) contains 15 objects that meet these criteria. 

2. Additionally, galaxies were drawn from the | O iii]+H/3 
narrowband excess catalog that we derived in Dressier 
et al. (2011b). The present star-forming sample contains 
1 1 objects that were selected in this way. 

3. Finally, we also detect R23 from ultra-faint emission lines 
that we originally considered Lya candidates, but, upon 
the follow-up observations described below, we determined 
that the discovery line was [Om] 74959, [Om] 75007, 
or H/3. Six galaxies with measurable R23 fall under this 
classification. 

In practice, the objects described in Class 3 are difficult to 
evaluate; their imaged counterparts are sometimes ambiguous 
due to their faintness and uncertain position within the blind- 
search slit. Therefore, in this paper we focus on the former two 
classes. For completeness, in Tables 2 and 3 we list the subset 
of these objects where R23 could be measured; however, we do 
not consider them further in this paper. 

The DEIMOS observations of these follow-up slit masks were 
carried out in 2011 January and 2012 January. A summary 
of observations for each follow-up mask is given in Table 1. 
In four of the five masks, we used slit widths and P.A.s that 
were matched to the search data (1"5 wide and 90° east of 
north). On mask L, in order to better locate objects detected 
through blind spectroscopy, we used a slit orientation that is 
orthogonal to the Venetian blind search slits. This method also 
allowed for narrower slits (1"2). All observations were carried 
out using the GG495 blocking filter with 830G grating. Under 
this configuration, we achieved a spectral resolution of 3.7 A 
(2.9 A) for a source that fills the 1"5 (1"2) slit. The grating angle 
was chosen to give a central wavelength of 7270 A, with blue 
coverage down to 5500 A for most slits. 

The DEIMOS spectra were reduced using the DEEP2 
DEIMOS data reduction pipeline (Cooper et al. 2012), with 
an updated optical model for the 830G grating (P. Capak 2010, 
private communication). The data were flux calibrated using 
observations of several spectrophotometric standard stars, taken 
through a 1'.'5 slit at the parallactic angle. The stars used were 
G191B2B, GD50, Feige 66, Feige 67, and Hz 44 (Massey 
& Gronwall 1990; Oke 1990), and the data for each were 
taken from the ESO spectrophotometric standard star database. 
Sensitivity functions derived from these stars differ systemati- 
cally between observations, with offsets up to 30%. After shift- 
ing the lower sensitivity functions to match the highest one, the 
observations agree at the 2%-3% level, indicating an excellent 
relative calibration. 

Finally, we have verified that the effects of differential 
atmospheric refraction have a negligible impact on our measured 
line flux ratios. In order to quantify possible differential slit 
losses, we calculate the component of the atmospheric refraction 
that falls perpendicular to the slit in each observation (frame) 
of each mask. The worst case occurs on mask L, where the slits 
are narrowest and the slit-P.A. was the furthest from parallactic. 
For these data, the refraction (perpendicular to the slit) between 
6000 A (near [On]) and 8000 A (near [Om] and H/3) ranges 
from (y/04 to 0H33 (Fillipenko 1982). Therefore, since the 


6 Because the COMSOS photometric redshifts include narrowband 
photometry and use templates that include emission lines, the photometric 
redshifts are often very close to the spectroscopic redshifts for the 
emission-line-selected galaxies in our survey. 
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Table 1 

DEIMOS Follow-up Observation Summary 


Date 

Mask Name 

Mask R.A. 
(J2000) 

Mask Deck 
(J2000) 

Mask P.A. 
(deg) 

Slit P.A. 
(deg) 

Slit Widths 
(") 

Exposure Time 
(hi) 

Seeing 

(") 

201 1 Jan 27 

D 

10:00:22.97 

02:09:28.9 

85 

90 

1.5 

6.5 

0.6 

201 1 Jan 28 

F 

10:01:11.46 

02:10:27.4 

106 

90 

1.5 

6.3 

0.8 

2012 Jan 22 

M 

10:00:24.25 

02:05:19.8 

85 

90 

1.5 

4.9 

1.0 

2012 Jan 22/24 

L 

10:00:22.56 

02:15:29.8 

13 

0 

1.2 

5.3 

0.9 

2012 Jan 23/24 

Q 

10:00:28.96 

02:20:05.1 

95 

90 

1.5 

6.8 

1.0 


Notes. Coordinates, position angles (P.A.s), exposure times, and seeing are given for each observed mask. On masks D, F, M, and Q, the slit P.A.s and 
slit widths were matched to the IMACS Venetian blind spectroscopic data (described in Dressier et al. 2011b and Henry et al. 2012). For mask F, the 
slits were rotated 90° relative to the search data. 


guiding was done in the R band (in between the observed 
wavelengths of the emission lines), we estimate that the slit 
losses on the red and blue parts of our spectra differ by no 
more than 4% (for 1" seeing). In most frames, the effect is even 
smaller, so we conclude that no systematic correction is needed 
to interpret our emission-line ratios. 

2.2. COSMOS Imaging Data 

In order to derive stellar masses (Section 3), we use the 
wide range of imaging data provided by the COSMOS team. 
A more detailed description of the catalogs can be found in 
Capak et al. (2007) and Ilbert et al. (2009). In summary, the 
photometric data include Galaxy Evolution Explorer ( GALEX ) 
imaging, a wealth of ground-based, broadband optical and near- 
infrared imaging, and the four Spitzer/IRAC bands. Further- 
more, our spectral energy distribution (SED) fits are improved 
by the inclusion of intermediate-band imaging and updated 
Subaru/SuprimeCam z'-band data (from observations made 
with the new, fully depleted, red-efficient CCDs). Since 
the optical and near-infrared photometry is performed in 
3" apertures on data with homogenized point-spread func- 
tions (PSFs) (1"5 FWHM), we apply a point-source aper- 
ture correction of —0.28 mag to all optical and near- 
infrared bands. For the IRAC data we use 3" 8 diameter 
apertures and apply point-source aperture corrections of —0.29, 
—0.33, —0.51, and —0.59 mag to bands one through four. No 
aperture corrections are applied to the GALEX data, as these 
magnitudes were derived from fits to the PSF. Finally, we note 
that the recommended zero point offsets have been applied 
(Capak et al. 2007; Ilbert et al. 2009), and a Galactic foreground 
extinction correction is made using the Cardelli et al. (1989) 
extinction curve with E(B — V) = 0.019. 

3. ANALYSIS 

3.1. Emission-line Measurements 

In total, we identify 47 galaxies at z ~ 0.6-0. 7; of these, 34 
have complete coverage of the R23 ([O m |+[0 u]/H/i; Pagel 
et al. 1979) metallicity indicator. Their spectra are shown in 
Figure 1. The remaining 13 galaxies at these redshifts either 
have essential lines lost in the OH airglow spectrum or are Class 
3 objects (described above), where the observed spectrum is 
especially faint because the galaxy may not be centered in the 
slit. Only six Class 3 objects (without mass measurements) have 
R23 measurements; their spectra are shown in gray at the bottom 
of Figure 1 . Excluding these objects leaves 28 galaxies with both 
mass and R23 constraints. In Section 3.3 we show that two of 
these galaxies may contain active nuclei, so we adopt a final 
sample that includes the remaining 26 star-forming objects. 



3720 3740 4860 4880 4960 5000 

Kest 


Figure 1. The spectra of 34 galaxies with complete coverage of the R23 
metallicity indicator ([Oiii]+[Oii]/H/?; Pagel et al. 1979) are shown. The 
spectra of 26 objects that comprise our MZ sample are shown in black and 
are sorted by stellar mass with higher masses at the top (the same as their 
appearance in Tables 2 and 3). The two AGN candidates are shown in gray at 
the top and also continue the mass-ordered scheme. The Class 3 objects with 
ambiguous optical counterparts are shown at the bottom (in gray) in the same 
order as their appearance in the tables. All spectra are normalized by their flux in 
[O hi] 7.5007, so it is apparent that higher-mass galaxies have relatively stronger 
H/3 and [O n] emission. The strong “absorption” trough near 4845 A in one 
object (980802) is a detector artifact. 

Emission-line fluxes are measured for all of the z ~ 0.6-0. 7 
galaxies that we identified. Gaussian profiles are fit to each 
line, and uncertainties are determined through a Monte Carlo 
simulation where noise is added to the spectrum and the 
fits are repeated. For [On], the two doublet components are 
fit simultaneously. As an example, Figure 2 shows a typical 
spectrum and the Gaussian fits to the emission lines. 
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Figure 2. The DEIMOS spectrum of object 54.5+4-0. 18 measures the R23 oxygen abundance indicator with high signal to noise. The Gaussian fits used to measure the 
emission-line fluxes are shown in red. From left to right, the lines shown are [O n] XX3121, 3729, H/3, [O hi] A4959, and [O hi] A.5007. This spectrum is representative 
of the fainter, lower-mass galaxies in our sample (see Tables 2 and 3). 

(A color version of this figure is available in the online journal.) 


To facilitate a correction for H/i stellar absorption, we 
also measure emission-line equivalent widths. Because it can 
be difficult to assess the uncertainty on the continuum flux 
density when it is measured spectroscopically, we compare 
two methods. On one hand, we measure the equivalent widths 
directly from the spectra. In 25/28 objects in our sample, we 
detect continuum, although it is sometimes weak. As a second 
method, we compare the emission-line fluxes from the DEIMOS 
data to the continuum under the emission line from the SED 
fit (described below). After accounting for a systematic offset 
(because the emission-line fluxes are subject to slit losses but 
the SED fits use total fluxes), the equivalent widths measured by 
these different techniques agree to within 50%. We adopt this 
level of uncertainty on the H fi equivalent widths. 

Finally, having measured the fluxes and equivalent widths 
of our emission lines, we apply a correction for the stellar 
absorption. Because the stellar absorption component is much 
broader than the emission line, this correction depends on the 
spectral resolution. At the resolution of our data, the correction 
is approximately 1 A in equivalent width (Cowie & Barger 2008; 
Zahid et al. 2011). While the amount of absorption also depends 
slightly on the age of the stellar population, the main source of 
uncertainty in this correction is from the observed H/i emission 
equivalent width, as discussed above. We propagate this error in 
our H/i fluxes and the metallicities that we infer in Section 4. 

3.2. Stellar Masses, SFRs, and Dust 
Constraints from SED Fitting 

Stellar masses are determined by fitting template stellar 
populations to the broad- and intermediate-band photometry 
described in Section 2. For this task, we use FAST (Fitting 
and Assessment of Synthetic Templates; Kriek et al. 2009). 
The choice of population synthesis templates can have an 
important impact on the properties derived from SEDs, because 
the contribution from the (infrared-bright) thermally pulsing 
asymptotic giant branch (TP-AGB) stars differs from model 
to model. The most widely used templates are Bruzual & 
Chariot (2003, BC03), Maraston (2005), and Chariot & Bruzual 
(2007; Bruzual 2007); the latter two options provide a larger 
contribution from TP-AGB stars (and correspondingly smaller 
masses). Nevertheless, the proper contribution from TP-AGB 


stars remains a subject of debate (Kriek et al. 2010; Conroy & 
Gunn 2010; Zibetti et al. 2013). Therefore, in order to facilitate 
comparisons with the literature, we derive stellar masses using 
the BC03 models. We have verified that the Chariot & Bruzual 
and Maraston et al. models produce stellar masses that are 
systematically smaller by approximately 0.1 dex. 

In order to assure accurate stellar population constraints, it is 
also important to account for emission-line contribution to our 
SEDs. (The stellar synthesis templates that we fit to the SEDs 
do not contain nebular features. Failing to remove emission- 
line contamination can result in incorrectly inferred ages and 
stellar masses; Schaerer & de Barros 2009; Atek et al. 2011.) 
We use our line fluxes to calculate the small contribution from 
[O ii] to the r band, and [O m] and H/i to the i band. For the 
galaxies in our sample, we find that l%-5% of the /--band light 
and 2%-15% of the /-band light can be attributed to emission 
lines. These contributions are subtracted from the r- and /-band 
continuum flux densities. Emission from Ha, on the other hand, 
falls between the z! and J bands, so it does not contaminate 
any of the photometry in the present analysis. Finally, we 
consider emission-line contamination to the intermediate-band 
photometry. In these cases, the corrections will be larger and 
more uncertain. Therefore, we exclude the band at 624 nm that 
includes [Oil] and the band at 827 nm that covers [Om] and 
H/i. The remaining 10 intermediate bands should be relatively 
unaffected. 

The grid of stellar population parameters that we fit with 
FAST includes: a set of exponentially declining star formation 
histories with e-folding times, r, ranging from 40 Myr to 10 Gyr; 
characteristic stellar population ages ranging from 50 Myr to the 
age of the universe at z ~ 0.6-0. 7; and Ay — 0-3 for a Calzetti 
et al. (2000) extinction curve. Additionally, we use a Chabrier 
(2003) IMF with metallicities of 0.004, 0.008, and 0.02 (solar). 
Supersolar metallicities are excluded, because, as we will show 
in Section 4, the galaxies in our sample mostly have subsolar 
to solar gas-phase metallicities. The derived stellar masses, star 
formation rates (SFRs), and visual extinction values are given in 
Table 3. Because of the intermediate-band photometry, the dust 
extinction constraints exclude much of the allowed parameter 
space. Therefore we use the SED-derived dust constraints to 
correct our emission-line ratios, including these uncertainties in 
our error budget. In making this correction, we also account for 
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Figure 3. The [O iii]/H/1 ratio as a function of stellar mass can be used to identify 
AGNs, similar to the [Oiii]/H/3 vs. [Nn]/H« "BPT” diagram (Baldwin et al. 
1981; Juneau et al. 2011). Red points show our emission-line-selected sample, 
while contours (with levels defined arbitrarily) show the SDSS data. The solid 
black lines show the demarcation between star-forming galaxies (below and 
to the left), objects with active nuclei (above and to the right), and composite 
objects (the region in between the black lines around 10 10 — 10 11 Mg). The two 
highest mass objects in our sample, 55.5-3-0.60 and 73.5+5-0.54, fall outside 
of the z = 0 star-forming locus, so we identify them as candidate AGNs. 

(A color version of this figure is available in the online journal.) 

the fact that nebular extinction is 2.3 times higher than stellar 
extinction on average (Calzetti et al. 2000; Cresci et al. 2012). 
While the precise relation between stellar and nebular extinction 
remains controversial (Cowie & Barger 2008; Cresci et al. 2012; 
Wofford et al. 2013), we adopt the “standard” Calzetti et al. 
(2000) relation to allow comparison with other studies. 

In Sections 6 and 7 we draw conclusions based on the SFRs 
that we have derived from these SED fits. In order to verify that 
the SFRs are reliable, we have carried out an assessment of the 
systematic uncertainties by comparing SFRs derived from [O n], 
H/3, and the SED fits. The results, outlined in the Appendix, 
show that the agreement between the different diagnostics is 
good; systematic offsets are smaller than 0.2 dex. This level of 
uncertainty does not affect our conclusions. 

3.3. Contamination from AGNs 

In order to measure the MZ relation of star-forming galaxies, 
it is important to ensure that the observed emission lines do 
not originate from an active galactic nucleus (AGN). The 
typical approach for low-redshift galaxies is to use the [Om]/ 
H/J and [Nn]/Ha emission line ratios (the canonical BPT 
diagram; Baldwin et al. 1981). However, for the redshift of our 
sample, Ha and [N ii] fall at infrared wavelengths. Therefore, 
in lieu of follow-up spectroscopy, we turn to an alternate 
diagnostic: the Mass-Excitation (MEx) diagram (Juneau et al. 
2011). This approach, shown in Figure 3, uses stellar mass 
as a proxy for the Ha /[Nil] A6583 ratio (relying on the MZ 
relation). The resulting diagram appears qualitatively similar 
to the traditional BPT diagram, with AGNs falling toward 
the top and right of the plot. In this figure, the data for 
our galaxies (red points) are compared to those from the 
SDSS (contours), and the solid black line demarcates the 
difference between star-forming galaxies, AGNs, and composite 
objects (as defined by Juneau et al. 201 1). While the evolution 
of the MZ relation compromises somewhat the use of this 


method, we note that two of our galaxies are found in the 
area identified as a possible AGN. Since metallicity evolution 
works in the sense of shifting the proxy relation between 
Ha/[N n] 16583 and mass, it shifts the star-forming locus (and 
black threshold curves) to the right in Figure 3. Therefore, the 
two objects in the AGN or composite part of the sample may in 
fact be part of the non-AGN distribution; however, we take the 
conservative step of excluding them from further analysis. We 
list the measured and derived properties of the two candidate 
AGNs at the ends of Tables 2 and 3. 

4. CALCULATING OXYGEN ABUNDANCES 

Interpreting metal abundance measurements requires that we 
account for systematic uncertainties. Different methods for mea- 
suring strong oxygen abundances yield results that are offset by 
up to 0.7 dex (Kewley & Ellison 2008; Lopez-Sanchez et al. 
2012; Andrews & Martini 2013). On one hand, photoioniza- 
tion models that are used to derive theoretical calibrations (i.e., 
Kewley & Dopita 2002) are often based on simplistic assump- 
tions about Hii region geometries and ionizing spectra (van 
Zee et al. 2006). On the other hand, empirical methods, which 
correlate line ratios with Hu region electron temperatures (i.e., 
Pettini & Pagel 2004; Pilyugin & Thuan 2005), generally give 
lower metallicities than photoionization model-based calibra- 
tions (McGaugh 1991; Kewley & Dopita 2002). Because elec- 
tron temperature measurements may be overestimated when 
temperature and density gradients are present in Hii regions, 
it is possible that the empirical calibrations are biased toward 
low metallicities (Peimbert & Costero 1969; Stasinska 2005). 
Alternatively, Nicholls et al. (2012) have suggested that discrep- 
ancies between electron-temperature measurements and theo- 
retical strong-line estimates can be explained if the electrons in 
H ii regions deviate from an equilibrium Maxwell-Boltzmann 
distribution. 

Ultimately, the differences between metallicity calibrations 
are still not fully understood, but the systematic offsets can be 
accounted for using a set of transformation equations given by 
Kewley & Ellison (2008). In the sections that follow (except for 
Section 6), we use the calibration from Kobulnicky & Kewley 
(2004, hereafter KK04), which is the mean of two theoretical 
R23 calibrations in the literature (McGaugh 1991; Kewley & 
Dopita 2002). In Section 6 we instead use the calibration from 
Maiolino et al. (2008) so that we may compare our data to results 
from Mannucci et al. (201 1). 

4.1. High Metallicity or Low? Determining the Branch ofR23 

Metallicities derived from R23 can be degenerate, since this 
diagnostic is double valued. At high metallicities, the oxygen 
lines are weaker because cooling is efficient and the H ii regions 
have lower electron temperatures. On the other hand, at low 
metallicities, the overall decrease in oxygen relative to hydrogen 
also imprints a decrease in R23 with decreasing metallicity. The 
“turnover” metallicity that demarcates the transition between 
the upper and lower branches of R23 depends on the ionization 
parameter and differs among the calibrations that are in the 
literature. For the typical ionization parameters of our galaxies 
and the KK04 calibration, the turnover metallicity is around 


7 The ionization parameter is defined as the ratio of the ionizing photon 
density to the hydrogen density. It can be written as U = Q/4nr 2 n j/c, where 
Q is the ionizing photon rate, r is the radius of the H n region, n is the 
hydrogen density, and c is the speed of light. The parameterization q = U x c 
is also commonly found in the literature. 
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Table 2 

Emission-line Measurements 


Henry et al. 


IMACS ID 

COSMOS ID 

R.A. 

(J2000) 

Deck 

(J2000) 

z 

F([0 hi] A5007) 

F([Om] A.4959) 

F(UP) 

F([On] A3727) 

EW(H/3) 
A (rest) 

27.5+5-0.47 

799604 

10:00:11.314 

+02:04:07.60 

0.675 

16.7 ±6.3 


39.2 ±0.7 

53.9 ± 1.3 

12.7 

99.5-3-0.84 

1235118 

10:00:53.842 

+02:22:06.88 

0.685 

19.75 ± 12.7 


25.1 ±0.6 

58.0 ± 1.2 

8.5 

90.5+6-0.63 

1265343 

10:00:09.101 

+02:19:52.02 

0.678 

10.2 ±2.7 


6.3 ±0.3 

16.4 ±0.4 

9.2 

100.5-2-0.41 

1234271 

10:00:50.551 

+02:22:22.66 

0.683 

13.4 ±2.2 


19.0 ±0.4 

48.5 ±0.9 

9.7 

65.5-1-0.69 

1003946 

10:00:47.710 

+02:13:36.83 

0.620 

26.0 ± 0.4 


19.4 ±3.0 

40.9 ±0.7 

9.6 

87.5+4-0.86 

1242162 

10:00:18.502 

+02:19:06.82 

0.656 

15.4 ±0.7 

5.4 ± 0.4 

5.4 ± 1.2 

19.4 ±0.8 

9.2 

95.5+4-0.94 

1236769 

10:00:21.331 

+02:21:05.22 

0.667 

158.5 ± 1.4 

52.5 ±0.5 

40.8 ±0.5 

112.3 ± 1.0 

32.0 

35.5+5-0.23 

794540 

10:00:09.038 

+02:06:08.16 

0.620 

57.0 ± 0.4 

19.1 ±0.4 

19.6 ±0.5 

70.5 ± 1.0 

14.8 


1015084 

10:00:41.863 

+02:09:12.66 

0.678 

309.7 ± 0.9 

104.1 ±0.5 

75.1 ±0.4 

200.3 ± 1.0 

39.5 

54.5+6-0.09 

1036285 

10:00:01.553 

+02:10:51.89 

0.684 

53.3 ± 1.0 

17.2 ±0.6 

18.5 ±0.3 

49.6 ±0.6 

27.0 

25.5+7-0.08 

801017 

09:59:58.039 

+02:03:38.89 

0.638 

36.9 ± 0.4 

12.8 ±0.4 

13.2 ±0.8 

46.9 ± 1.0 

12.1 

23.5+2-0.33 

775548 

10:00:32.287 

+02:03:09.25 

0.670 

32.8 ± 0.7 

12.2 ± 1.3 

14.8 ±0.5 

44.8 ± 1.1 

22.3 


1035959 

09:59:59.347 

+02:10:57.10 

0.628 

35.2 ±0.3 

12.7 ±0.2 

13.2 ±0.7 

34.7 ±0.6 

11.7 


982000 

10:01:07.140 

+02:12:08.99 

0.621 

235.1 ±0.6 

77.2 ± 0.5 

53.3 ±0.7 

112.7 ±0.9 

44.0 


1013478 

10:00:18.701 

+02:09:41.74 

0.635 

34.9 ± 0.3 

12.1 ±0.3 

11.2 ±0.5 

40.7 ± 0.6 

28.4 


980802 

10:01:07.135 

+02:12:38.22 

0.621 

27.3 ± 0.4 

9.2 ±0.6 

11.9 ±0.7 

40.5 ±0.7 

49.3 

28.0+7. l a 

799190 

10:00:03.035 

+02:04:13.98 

0.638 

167.5 ± 0.6 

57.3 ±0.4 

35.2 ±0.8 

91.3 ± 1.0 

56.3 


770439 

10:00:42.334 

+02:05:43.53 

0.629 

34.4 ± 0.4 

11.5 ±0.4 

9.8 ±0.7 

23.9 ±0.8 

22.9 


1009808 

10:00:44.558 

+02:11:18.04 

0.616 

96.9 ± 0.4 

32.5 ± 0.5 

22.2 ±0.9 

45.2 ±0.8 

38.0 


1237667 

10:00:52.368 

+02:20:58.18 

0.640 

59.2 ± 0.4 

21.2 ±0.4 

13.6 ±0.4 

24.9 ± 1.0 

40.9 


772773 

10:00:25.567 

+02:04:44.36 

0.616 

21.5 ±0.5 

8.2 ±0.8 

6.9 ± 1.0 

22.3 ± 1.0 

11.3 

48.5-4-0.66 

989316 

10:01:08.537 

+02:09:21.70 

0.674 

13.8 ±0.5 

4.7 ± 1.0 

4.96 ±0.2 

16.1 ±0.7 

17.7 

54.5+4-0.18 

1010994 

10:00:16.097 

+02:10:51.60 

0.634 

28.7 ± 0.2 

10.0 ±0.3 

5.4 ±0.5 

13.4 ±0.5 

12.2 

47.5+3-0.85 

1015516 

10:00:27.792 

+02:09:06.43 

0.679 

34.1 ±0.3 

12.4 ±0.2 

8.5 ±0.2 

17.7 ±0.6 

43.5 


1010583 

10:00:21.408 

+02:11:00.43 

0.639 

45.2 ± 1.5 

14.6 ±0.8 

10.6 ±2.1 

28.5 ±0.6 

27.9 


1217842 

10:00:54.658 

+02:19:01.81 

0.633 

21.4 ±0.3 

8.4 ± 0.4 

5.1 ± 1.7 

9.8 ±0.8 

19.6 

Objects lacking clear optical counterparts in the COSMOS data 

52.5+8-0.69 


09:59:51.374 

+02:10:23.28 

0.678 

15.6 ±0.8 

3.7 ±0.9 

2.1 ±0.1 

3.6 ±0.4 


45.5+7-0.25 


09:59:55.699 

+02:08:38.51 

0.678 


2.2 ± 0.7 

2.6 ±0.4 

9.6 ± 1.0 

6.5 

46.5-7-0.04 


10:01:25.152 

+02:08:52.34 

0.672 

6.2 ±0.5 

2.5 ± 1.5 

2.4 ±0.8 

5.4 ±0.7 


65.5-2-0.14 


10:00:50.904 

+02:13:36.52 

0.624 

14.1 ±0.3 

5.0 ±0.9 

5.4 ±0.2 

6.8 ±0.8 


47.5+6-0.90 


10:00:07.145 

+02:09:07.81 

0.639 

14.4 ±0.4 

5.0 ±3.8 

3.4 ±0.8 

<9.5 


87.5-1-0.06 


10:00:48.792 

+02:19:06.54 

0.620 

15.4 ±0.3 

5.3 ± 1.6 

6.0 ± 1.7 

1 .9 ± 1.3 


Candidate AGN 

55.5-3-0.60 

984577 

10:01:01.126 

+02:11:07.85 

0.623 

39.5 ± 0.4 

13.7 ±0.4 

16.4 ±0.5 

51.6 ±0.7 

9.9 

73.5+3-0.54 

998509 

10:00:25.046 

+02:15:38.15 

0.642 

38.7 ±0.5 

12.4 ±0.6 

48.4 ±0.8 

93.3 ±0.9 

9.4 


Notes. Emission-line fluxes are in units of 10 18 erg s’ 1 cm 2 . Equivalent widths and HfS fluxes are given as measured (without applying the 1 A stellar absorption 
correction). Objects lacking an IMACS ID were selected strictly from the COSMOS narrowband (NB816) imaging, as we described in Section 2. The [O ii] flux is the 
sum of the AA3726, 3729 doublet fluxes. In six cases, one of the lines in the [O hi] doublet lines fell on an OH sky line; in this situation, metallicities are derived by 
assuming a flux ratio of f 1 ([Om] A5007)/F(|0 ill ] A4959) = 3. 

a This galaxy was drawn from our earlier (wider and shallower) survey first presented in Martin et al. (2008). 


12 4- log(0/H) ~ 8.4, which is reached when log(R23) ~ 
0.9-1. 0. 

Many different methods have been suggested for breaking 
the degeneracy between the upper and lower branch. The most 
rigorous of these is to compare to an alternative metallicity 
diagnostic. Kewley & Ellison (2008) advocate for the use of the 
[N ii] A6583/[Oii] or [N ii] A6583/Ha ratios to first provide a 
metallicity estimate that identifies the branch. However, this 
method is more challenging at z > 0.5 (requiring infrared 
spectroscopy to reach very faint [N n] AA6548, 6583 lines) and 
impossible from the ground at z > 2. The R23 diagnostic, 
on the other hand, can be measured from the ground out to 
z ~ 3 (Maiolino et al. 2008). It is therefore imperative that 
we learn how to break the R23 -metallicity degeneracy without 
observations of the [N ii] 7.7.6548, 6583 lines. Here, we examine 
other methods for establishing the R23 branch that have been 
proposed in the literature. 


In Figure 4 we investigate the branch of R23 by plotting 
four different diagnostics as a function of stellar mass. For 
comparison, we show the four H/i -detected galaxies (0.60 < 
z < 0.85) from the PEARS Survey (Probing Evolution and 
Reionization Spectroscopically; Xia et al. 2012). In the upper 
left panel, we show that R23 decreases with increasing stellar 
mass for most of our sample. Furthermore, three of the four 
galaxies from Xia et al. follow the trends from our emission- 
line-selected sample. This diagnostic paints a clear picture: if 
metallicity decreases as a function of decreasing mass, most of 
the galaxies in our sample must fall on the upper branch of the 
R23 indicator. The two lowest mass galaxies shown in Figure 4 
(1217842 from our sample, and 246 from Xia et al.) hint at a 
possible turnover at M ~ 8. 0-8. 5, indicating that these galaxies 
may fall on the lower branch of R23. 

Because it is not ideal to introduce a mass dependence 
on our metallicity derivations, we explore other emission-line 
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Table 3 

Derived Properties 


Henry et al. 


IMACS ID 

COSMOS ID 

M b 

U-B 

log M t 

log SFR 

Ay 

log R23 

log 032 

12+log(0/H) 

log? 

27.5+5-0.47 

799604 

-20.93 

0.52 

10 05 +0 09 
iu.uj_o 03 

1 70 + 0.34 
A ‘ ZJ — 0.21 

1-80-o.4 3 8 

0.68 + 0.10 

-0.94 + 0.16 

8.80/° 0 08 

7 00 + ° 09 
7-UU -0 .13 

99.5-3-0.84 

1235118 

-20.71 

0.45 

n 77 + 0.04 
y - 1 '- 0.00 

n 95 + 0.20 
- 0.00 

l-40 + _° 0 ° 1 | 

0.77 ± 0.07 

-0.78 ± 0.22 

8.63/g 0 ) 8 

7 05 +0 - 14 
7 - u - 0.10 

90.5+6-0.63 

1265343 

-19.24 

0.61 

Q 74+O.OO 
y ' 7 —0.25 

—0 S 4 +012 
u - Ji t- 0.00 

0 20 +0 30 
U-ZU — 0.05 

0.63 ± 0.08 

-0.14 + 0.11 

Q QQ+0.07 

o.oo_0.06 

7 65 +0,12 
7 ■°' 5 — 0.07 

100.5-2-0.41 

1234271 

-20.39 

0.42 

Q A7+0.08 

^•oz _ 009 

0 62 +0-38 
U-OZ — 0.05 

1 oo +0 - 28 

i.UU— 0 3i 

0.71+0.10 

-0.75 + 0.19 

8-76^0 

711 +0.06 
7 - A 1 - 0.10 

65.5-1-0.69 

1003946 

-19.56 

0.42 

q 17+O.O9 
yL '-0.07 

—0 4S+0- 44 
0.06 

n An+O-O^ 
°- 6 °- 0.21 

0.62 ± 0.09 

-0.26 ± 0.05 


7 ^^+0.09 
7 - JJ - 0.11 

87.5+4-0.86 

1242162 

-19.15 

0.44 

9 15 -006 

—0 42 +0 07 
u -^ z _o.oo 

0 oo +018 

u -uu_ 000 

0.78 + 0.12 

0.03 ± 0.04 

8 . 68 /q 18 0 

7 66 +0-14 
/.oo_o 12 

95.5+4-0.94 

1236769 

-20.08 

0.37 

n 1 9 + 0.08 
y ‘ AZ — 0.05 

0 30 +0 ' 59 
u.JU -0 22 

0 60 +0 - 27 
— 0.40 

0.85 ± 0.05 

0.09 + 0.10 

O 41+O.I6 

°-^ A -o.io 

7 C 7 + 0 . 14 
7 - J 7 - 0.12 

35.5+5-0.23 

794540 

-19.89 

0.30 

q A9+O.04 

^• uz -0.04 

0 02 +0,22 
uuz - 0.05 

0 2O +0 - 40 
u ' zu -0.07 

0.85 ± 0.06 

-0.03 ± 0.07 

0 co+0.09 

-0.06 

7 56 +0 - 07 
/ - JO - 0.12 


1015084 

-20.26 

0.21 

8-99!°o 0 0 9 o 

0-33!° 0 ' 1 o 5 0 

0-001ZI 

0.89 + 0.01 

0.32 ± 0.02 

0 ^7+0.01 
°' J/ -0.02 

7 00 + 0.02 

/• 0-3 _0.03 

54.5+6-0.09 

1036285 

-18.93 

0.40 

8-98$£ 

0.08!° 0 3 4 8 6 

0-60/° 0 4 6 4 o 

0.87 ± 0.08 

-0.03 + 0.14 

0 ^7+0.02 

7 54 +0 - 14 
7 ■ J ^-0.05 

25.5+7-0.08 

801017 

-19.31 

0.37 

8-9+°o°03 

—0 06 +0 0 ° 
U.UO_o25 

0 20 +0 - 22 
U ' ZU -0.20 

0.83 ± 0.06 

-0.04 ± 0.06 

0 f . 1+0.10 
O-Ol_ 008 

7 56 +0 09 
7 - JU -0.08 

23.5+2-0.33 

775548 

-19.32 

0.34 

8 . 86 Z ° 0 ° 0 6 3 

—0 7^+018 
-7-0.16 

n 9o +014 
u>zu -0.20 

0.78 ± 0.04 

-0.06 ± 0.05 

0 77+0.05 
7 2 — 0.04 

7 61 +0 09 
/•O 1-0.08 


1035959 

-19.11 

0.37 

0 0 'j+0.09 
-0.09 

—0 36 +0 - 20 

u --70_o is 

0 20 +0 - 23 
u ' zu -0.20 

0.75 ± 0.06 

0.08 ± 0.06 

0 7C+0.07 
7 —0.06 

7.74/°/° 


982000 

-19.54 

0.25 

8.82/q° 1 5 4 

0.43 !° 0 2 4 3 

0 40+0.24 
u -^ u -0.26 

0.92 ± 0.03 

0.32 + 0.08 

0 c 1 +0.08 

0.04 

i-so+Zl 


1013478 

-18.64 

0.36 

O 77 +O.O 8 
7 —0.08 

-0-28^ 

0 6 O +0 ’ 46 
u.ou _0 24 

0.96 ± 0.08 

- 0.12 + 0.11 

0 oc+0.21 

-o.i 6 

7 40+ 0 09 
7 -^ u -0.05 


980802 

-18.74 

0.37 

8-63i° 0 ° 0 9 8 

—0 56 +0-33 
-0.50 

0 20 +0 - 27 
u>zu -0.20 

0.83 ± 0.05 

-0.11 ±0.07 

0 /cc+0.10 
5.OD-0 07 

7 59 + 0.11 

28.0+7.1 

799190 

-19.07 

0.24 

8.62^“ 

-0.04^ 

0 20 +0 - 43 
U-Z — 0.07 

0.96 ± 0.03 

0.33 ± 0.07 

0 49 + 0.07 
O-^J- 0.10 

7 77+O.O7 
' 7 -0.13 


770439 

-18.46 

0.36 

0 ^'j+O.lO 

°- JJ -o.oo 

-1-08/° 0 4 o 6 o 

0 20 +011 
u - zu -0.20 

0.84 ± 0.04 

0.22 ± 0.04 

0 ^.a+0.07 
o-O4_ 008 

1 on+0.06 
/•oU_o 07 


1009808 

-18.62 

0.28 

0 45 + 0 . is 

—0 13 +018 
u> 1 — 0.35 

0 20 +0-27 
u>zu -0.20 

0.89 ± 0.03 

0.39 ± 0.07 

q ^7+0.07 
O-J/-0.04 

7 -9 l/ggly 


1237667 

-17.99 

0.30 

8-43^ 

—0 51 +0-76 
0.02 

0 oo +0 - 82 
U -W_ 0.00 

0.87 ± 0.04 

0.51 + 0.08 

8-63To6 

8.05 + _° 0 % 


772773 

-18.72 

0.28 

8-40$& 

—0 61 +0-30 
u - oa - 0.02 

0 oo +015 
u - uu -o.oo 

0.81+0.08 

0.12 + 0.04 

s.eeTVoe 

7 77+O.O8 
7 ‘ 7 Z — 0.06 

48.5-4-0.66 

989316 

-17.94 

0.35 

0 72Q+0.23 
QDy - 0.13 

—0 21 +0 - 37 
U-ZA — 0.52 

0 60 +0 - 40 
u - ou -0.60 

0.89 ± 0.09 

-0.13 + 0.14 

0 49 +O.O 2 

7 44+0.03 
7 -^-0.18 

54.5+4-0.18 

1010994 

-18.18 

0.40 

q tjq+0.07 

o-JO- 0.00 

_i 41 + 0.06 
A '^ A - 0.01 

0 oo +004 
u - uu -o.oo 

0.92 ± 0.06 

0.46 ± 0.02 

8.49 /° 0 ° 0 8 

7 9 1 +0.05 
7 — 0.05 

47.5+3-0.85 

1015516 

-17.82 

0.36 

Q O 7 +O .73 
—0.07 

0 - 46 -°i 3 i 0 i 

1 AH+0-30 
1 ' 60 — 1.54 

1.06 + 0.12 

-0.08 ± 0.25 

8.40/° 0 21 0 

7 36 +0 ‘ 22 
/ - 3O -0.19 


1010583 

-18.00 

0.26 

0 77+O.2O 
°- z2 -0.16 

-0.06/° 0 2 s 8 6 

0 60 +0 38 
UOU -0.56 

0.96 + 0.11 

0.13 + 0.13 

q oq+0.28 

5. JO -0 26 

7.6O/°0 1 1 5 7 


1217842 

-18.00 

0.20 

7 og+0-09 
7 yo - 0.06 

—0.48/^32 

0 oo +0-22 
u - uu -o.oo 

0.85 + 0.16 

0.49 ± 0.05 

8-40/q 4 q 

8.O5/°0 22 

Objects with ambiguous counterparts in the COSMOS data 

52.5+8-0.69 







1.03 + 0.04 

0.73 ± 0.05 



45.5+7-0.25 







0.74 + 0.12 

-0.04 + 0.13 



46.5-7-0.04 







0.77 + 0.16 

0.21 + 0.10 



65.5-2-0.14 







0.68 ± 0.03 

0.45 ± 0.06 



47.5+6-0.90 







0.86 + 0.14 

0.61 + 0.35 



87.5-1-0.06 







0.57 + 0.13 

1.04 + 0.29 



Candidate AGN 

55.5-3-0.60 

984577 

-19.78 

0.62 

10 2 < S +0 01 
1U.ZO_ 00 0 

0 A7+ 0 00 
U - O/ -0.16 

1 ?o +01 ° 
1 -ZU — 0.20 





73.5+3-0.54 

998509 

-21.23 

0.54 

1 0 23 +0 00 
1U.ZJ_ 0 _06 

0 62 +0,55 

1 4O +0 - 20 
1 -^ U _ 0.12 






Notes. These quantities are derived from COSMOS imaging and our measured line ratios. Magnitudes and colors are in AB units and are derived by integrating 
the best-fitting SED under the appropriate bandpasses. Stellar masses and SFRs are given in solar units and were derived using FAST (Kriek et al. 2009). For this 
derivation we used a Chabrier (2003) IMF and a Calzetti et al. (2000) extinction curve. The line ratios, metallicity, and ionization parameter ( q ) are derived, taking 
into account a 1 A (equivalent width) correction to the H/3 fluxes, as well as the dust correction and its associated uncertainty. Metallicities and ionization parameters 
are calculated using the Kobulnicky & Kewley (2004) calibration. A dust correction has not been applied to the objects that have ambiguous optical counterparts in 
the COSMOS imaging survey. 


diagnostics in Figure 4. First, we compare to the emission- 
line ratios of [Om]/H/J and log([Om]/[Oii]) = 032 as a 
function of stellar mass. These line ratios have been proposed 
to break the R23 degeneracy by Maiolino et al. (2008), with 
low-metallicity solutions being preferred when [ 0 1 1 1 ] / H /3 > 5 
and 032 > 0.5. Indeed, Figure 4 shows that the galaxies above 
these thresholds also prefer to have higher values of R23 (and 
lower metallicities on the upper branch). Nevertheless, we still 
see very little evidence for a substantial population of galaxies 
extending to metallicities below the R23 turnaround. The only 
object that passes both the 032 and [Oiii]/H/3 thresholds is 
object 246 from Xia et al. (2012). 


The inference of upper branch metallicities for 3/4 of the 
objects presented by Xia et al. (2012) contrasts with their lower 
branch assumption. These authors argued for lower metallicities 
because of the relatively high values of [O m]/HyS and 032 
for most of their sample. However, comparison to the R23 
versus mass plot (upper left) shows that for at least three of 
these galaxies, higher-metallicity solutions are more plausible. 
(Although the galaxies presented by Xia et al. 2012 may prefer 
higher 032 and ionization parameters than the remainder of 
our sample, metallicity is not a strong function of ionization 
parameter on the upper branch, so their galaxies still follow our 
R23-stellar mass correlation.) 
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Figure 4. Our measured values of R23 decrease with increasing stellar mass, implying that most of the galaxies in our sample fall on the upper branch of R23. 
This diagram is contrasted with three other diagnostics that have previously been used to identify whether galaxies have low or high metallicities (Maiolino et al. 
2008; Hu et al. 2009; Xia et al. 2012). These include the [O hi] XX4959, 5007/H/3 flux ratio, the rest-frame equivalent width of H/l, and the ratio 032 = log ([O hi] 
XX4959, 5007/[0n] A.A.3726, 3729). With an H/l EW of 352 A (rest), object 246 from Xia et al. is not shown in the EW(H/?) panel. 


Ultimately, caution is required when using [Oiii]/H/3 and 
032 to determine the branch of R23, because these quantities 
are dependent on the ionization parameter. Observations of high- 
redshift galaxies show that their line ratios are offset to high 
values of |Om]/H/J at a fixed [Nn]/Ha (Shapley et al. 2005; 
Erb et al. 2006; Hainline et al. 2009). This offset is usually 
interpreted as evidence for a systematically high ionization 
parameter (Brinchmann et al. 2008), suggesting that some 
galaxies deviate from the local metallicity versus ionization 
parameter correlation. In summary, attempting to break the 
R23 degeneracy using quantities that depend on the ionization 
parameter could erroneously indicate lower branch solutions. 

In Figure 4, we also investigate the use of H/i equivalent 
width to select low-metallicity galaxies. Kakazu et al. (2007) 
and Hu et al. (2009) have found that, for their emission-line- 
selected sample, galaxies are likely to have very low metallicities 
(indicated by detectable [O in] 14363 emission) when their H/3 
equivalent widths are more than 30 A (rest). However, 8/26 
galaxies in our sample meet this criterion, and we do not detect 
[O in] 14363 in any objects. What is more apparent in Figure 4 
is that a cut in H/3 equivalent width merely selects against 
galaxies with M > 10 90 M 0 . As we show in Section 5, the 
Hu et al. sample has, on average, higher H/i equivalent widths 
and somewhat lower luminosities than the present sample of 
emission-line objects. We infer that it is not straightforward to 
use H/i equivalent width to discriminate between upper and 
lower branch metallicities. 


Guided by Figure 4 we conclude that, in the absence of other 
metallicity diagnostics, the correlation between stellar mass and 
R23 is the preferred method to break the R23 degeneracy. This 
exercise shows that for z ~ 0.6-0. 7 galaxies, the maximum 
value of R23 is reached between 10 8 0 and 10 8 5 M e . Therefore, 
we adopt upper branch solutions above M — 10 8 ' 2 M 0 . The 
two lowest mass galaxies (1217842 from our sample and 246 
from Xia et al. 2012) are less certain. For 1217842 we adopt 
a metallicity in the turnaround region (12+log(0/H) = 8.4), 
with error bars denoting the upper bound of the high-metallicity 
solution and the lower bound of the low-metallicity solution. The 
lowest-mass object from Xia et al. (#246) is tentatively assigned 
to the lower branch because of its more extreme line ratios, 
although we caution that the H/J measurement could be in error 
because it is blended with [O in] in their low-resolution data. 
Follow-up spectroscopy can better constrain the metallicities of 
galaxies near the R23 turnaround by providing the [Nil] /Ha 
and [Nii]/[Oii] ratios (Kewley & Ellison 2008). At present, 
however, the conclusions drawn in the remainder of this work 
do not depend on these two galaxies because their metallicity 
errors are large. Oxygen abundances for the entire sample are 
listed in Table 3. 

Our inference of upper branch metallicities at M > 
10 8 2 -10 8 5 Mq is inconsistent with results reported by Zahid 
et al. (2011). In contrast to Figure 4, Zahid et al. find that the 
mean value of R23 turns over around M ~ 10 9 2 M 0 for their 
Z ~ 0.8 sample. They interpret the turnover as evidence for 
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Figure 5. The addition of our 26 low-mass galaxies constrains the low-mass portion of the MZ relation at intermediate redshifts. For comparison, we have re-calculated 
the metallicities presented in Xia et al. (2012), assuming — contrary to these authors — that the three highest-mass objects have metallicities on the upper branch. Other 
studies are shown for comparison and are (when necessary) converted to a Chabrier (2003) IMF and a Kobulnicky & Kewley (2004) metallicity. The dashed and dotted 
gray lines that follow the Tremonti et al. (2004) MZ relation show the 68% and 95% contours for the SDSS data. For reference, solar metallicity is 12+log (O/H) = 
8.69 (Allende Prieto et al. 2001). 

(A color version of this figure is available in the online journal.) 


lower branch metallicities. However, as can be seen in their 
Figure 4, the data are increasingly noisy below M < 10 9 ' 5 M q , 
and objects with R23 > 10 have been removed as candidate 
AGNs. The coupling of these effects may bias the lowest-mass 
bin toward lower mean values of R23 and mimic the effect of a 
turnover in the R23 versus mass correlation. We conclude that 
on average, galaxies at z ~ 0.6-0. 7 fall on the upper branch of 
R23 for M > 10 s 2 — 10 s 5 Mq. At higher redshifts, where metal- 
licity evolution is significant (i.e., Erb et al. 2006), we expect 
that the transition from the upper to lower branch of R23 will 
occur at higher stellar masses. 

5. RESULTS 

5.1. The Mass-Metallicity Relation 

Figure 5 presents the MZ relation for our 26 star-forming 
galaxies and compares it to both the local relation (Tremonti 
et al. 2004) and other intermediate-redshift measurements from 
the literature. In each case, we have taken care to (when 
necessary) convert stellar masses to a Chabrier (2003) IMF 
(using multiplicative constants given in Savaglio et al. 2005; 
Cowie & Barger 2008) and convert metallicities to a KK04 
calibration (using equations given in Kewley & Ellison 2008). 
In the mass range 8.5 ^ M 0 5) 9.0, our data give a median 
value of 12+log(0/H) = 8.63, with an rms scatter of 0.12 dex. 

The intermediate-redshift MZ relation has been measured for 
higher-mass galaxies by several authors (Savaglio et al. 2005; 
Cowie & Barger 2008; Lamareille et al. 2009; Zahid et al. 2011). 
Among the higher-mass portion of our sample (M > 10 9 0 Mq, 
where we can compare directly to previous results) there is 
good agreement with all of the published MZ relations. While 
the addition of our data does not distinguish between these 
relations, we note that the Zahid et al. sample stands out as the 
largest (1350 galaxies, as opposed to fewer than 100 galaxies 


in each of the others). These authors discuss the differences 
between Savaglio et al. (2005), Cowie & Barger (2008), and 
Lamareille et al. (2009). They conclude that sample selection 
effects and linear fits that are biased by outliers may account 
for the differences among these works. While these effects 
are likely important, we also note that different dust correction 
methods maybe significant. On one hand, Zahid etal. (2011) and 
Lamareille et al. (2009) calculate R23 from equivalent widths, as 
they are less sensitive to dust extinction. Although this method is 
subject to systematic uncertainties correlated with galaxy colors, 
Liang et al. (2007) show that this effect has a relatively weak 
impact on metallicity (—0.2 to 0.1 dex). Cowie & Barger, on 
the other hand, take the approach of determining dust extinction 
from SED fits. However, they do not assume a higher dust 
extinction for nebular compared to stellar light, as their lower- 
redshift sample suggests that these quantities are consistent. 
Comparatively, the dust correction adopted by Savaglio et al. 
(2005) is cruder; they use the same extinction (Ay = 2.1) 
for all galaxies. If this correction is accurate for the median 
of their sample, then at high (low) masses the dust correction 
will be underestimated (overestimated), and metallicities will 
be overestimated (underestimated). Qualitatively, this effect 
could explain their steeper MZ slope. Additionally, it is worth 
noting that our data are inconsistent with an extrapolation of 
the Savaglio et al. relation. In the sections that follow, we take 
the results from Zahid et al. (2011) as the best measurement 
of the high-mass, intermediate-redshift MZ relation. 

The addition of our low-mass data allows for new constraints 
on the evolution of the MZ relation. In Figure 5 we compare 
to the local relation measured from the SDSS (Tremonti et al. 
2004). In the mass range 10 8 5 < M/M Q ^ 10 9 °, we find a 
mean metallicity that is 0.12 dex lower than the local relation 
(at M = 10 s 75 Mq). This trend confirms that metallicity 
evolution is relatively slow from intermediate to low redshifts. 
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Figure 6. The luminosity-metallicity relation for the present emission-line- 
selected sample is compared to correlations from the literature. Notably, we 
compare to the measurements of the ultra-strong emission-line galaxies reported 
by Hu et al. (2009). Because these galaxies have detected [O hi] X4363 emission, 
we assume that they fall on the lower branch of R23. Their metallicities are 
recalculated using the KK04 calibration (assuming no dust, as justified by Hu 
et al.). 

(A color version of this figure is available in the online journal.) 


Additionally, as our data are consistent with an extrapolation of 
the results from Zahid et al. (201 1), we agree with their finding 
that metallicity evolution is more significant at lower masses 
than at higher masses. This trend is qualitatively consistent 
with downsizing in the later phases (z ~ 1-2 to today) of 
galaxy evolution. After high-mass galaxies have assembled most 
of their stars, they have also made most of the metals that 
they will ever produce. At the same time, these higher-mass 
galaxies are less effective at removing metals through winds, 
and they are also less efficient at accreting gas (Dekel et al. 
2009). As a result, higher-mass galaxies are less capable of 
changing their metallicities after they have assembled most 
of their stars. Under this scenario, we would expect lower- 
mass galaxies to evolve onto the present MZ relation at later 
cosmic times. Hence, extending observations to lower masses 
at all redshifts will provide essential insights into the physics of 
galaxy assembly and metal production. We compare our data to 
theoretical predictions for the MZ relation in Section 7. 

5.2. The Luminosity— Metallicity Relation 

In Figure 6 we show the metallicity-luminosity relation 
derived from our data. While the MZ relation may represent 
a more fundamental correlation than the metallicity-luminosity 
relation, the latter allows comparison to studies where stellar- 
mass measurements have not been possible. Here, we focus 
on the ultra-strong emission line (USEL) sample presented by 
Kakazu et al. (2007) and Hu et al. (2009). Like the galaxies 
in our sample, the USELs lie at intermediate redshifts and are 
emission line selected. However, the majority of the USELs are 
low-metallicity galaxies with detected [Oiii] 74 3 63 emission. 
In Figure 6, we compare the metallicity-luminosity relation of 
the USELs to our sample. In order to make a fair comparison, 
we have calculated the metallicities of the USELs using R23 
and KK04, assuming lower-branch solutions. No correction for 
dust is applied, which Hu et al. (2009) argue is appropriate for 
these faint, low-mass objects. It is worth noting that, under these 


assumptions, the majority of USELs have 12+ log(0/H) > 8.0 
and would not be classified as extremely metal-poor galaxies 
(12+log(0/H) < 7.65; Kniazev et al. 2003). However, this 
classification changes when direct method metallicities are used: 
in this case, the oxygen abundances reported by Hu et al. (2009) 
are 0.03-0.74 dex lower (0.4 dex on average). 

Figure 6 shows clearly that the USELs have lower metal- 
licities than the emission-line-selected galaxies in the present 
work. In fact, this trend is not surprising. Because of the dif- 
fering selections for the USELs and the present sample, the 
former have H/l equivalent widths that are two times larger 
(on average) than the latter. This difference translates to higher 
specific (mass-normalized) SFRs for the USELs. One possible 
explanation is that more actively star-forming galaxies will have 
/1-band luminosities that are high for their stellar mass. In this 
scenario, galaxies that fall along the MZ relation will move to 
the right when plotted in the luminosity-metallicity relation, as 
in Figure 6. However, another explanation is possible. Within 
the local MZ relation, galaxies with higher specific SFRs have 
lower metallicities and follow a mass-metallicity-SFR plane 
called the Fundamental Metallicity Relation (FMR; Lara-Lopez 
et al. 2010; Mannucci et al. 2010, 2011; Yates et al. 2012; An- 
drews & Martini 2013). Hence, the differing metallicities be- 
tween the USELs and the present sample support the existence 
of the FMR in the relatively unexplored low-mass, intermediate- 
redshift regime. We next provide the first quantitative test of the 
FMR for intermediate-redshift galaxies below 10 9 M Q . 

6. THE FUNDAMENTAL METALLICITY RELATION 

In the previous section, we showed that, for a fixed luminosity 
(and possibly mass), galaxies with higher equivalent width 
emission lines (i.e., the USELs) have lower metallicities. This 
result supports the idea that star formation drives some of the 
scatter in the MZ relation, even out to intermediate redshifts. 
For local galaxies at fixed stellar mass, objects with higher SFRs 
have lower metallicities than lower-SFR galaxies (Ellison et al. 
2008; Mannucci et al. 2010, 2011; Yates et al. 2012; Andrews & 
Martini 2013). Mannucci et al. (2010) describe this relation as 
a plane and refer to it as the FMR. Remarkably, Mannucci et al. 
report that high-redshift galaxies also lie on the plane (albeit 
at lower masses and higher SFRs than most low-z galaxies), 
implying that the physics governing the evolution of metal 
enrichment has not changed over cosmic time. This finding 
has recently been confirmed for a large sample of intermediate- 
redshift (z ~ 0.6) galaxies from zCOSMOS (Cresci et al. 2012). 
Nevertheless, tests for evolution of the FMR have not yet been 
extended below M ~ 10 9 M e . Here, we provide the first such 
test. 

In light of the significant differences between metallicities 
derived from various calibrations (Kewley & Ellison 2008), it is 
essential that we use metallicities that are derived consistently 
with those in Mannucci et al. (2010, 2011; the latter work 
extends the initial analysis to lower masses). These metallicities 
are based on the calibration from Maiolino et al. (2008), which 
is semi-empirical. The high metallicities are derived from the 
photoionization models of Kewley & Dopita (2002), whereas the 
lower metallicities are constrained from direct T e measurements 
reported by Nagao et al. (2006). While the Maiolino et al. 
calibration of R23 does not include a dependence on the 
ionization parameter, in other calibrations (i.e., KK04) this 
dependence is weak on the upper branch of R23, where the 
present sample lies. Therefore, for this section, we re-calculate 
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log M/M 0 

Figure 7. The fundamental metallicity relation, derived using the Maiolino 
et al. (2008) metallicity calibration, shows that low-mass galaxies with high 
SFRs have systematically lower metallicities. The galaxy 23.5+2-0.33, which 
strongly influences the weighted mean of the residuals (see text), is shown as an 
open green symbol. Top: the MZ relation is plotted, with galaxies color-coded 
in three bins of SFR. The curves show the local FMR reported by Mannucci 
et al. (2011). Bottom: residuals from the local FMR (data minus relation given 
in Equation (1)) show increased scatter toward low stellar masses, owing to 
larger measurement errors for faint galaxies. (Note that residuals are calculated 
using measured SFRs rather than the mean values where Equation ( 1 ) has been 
evaluated for the top panel. As a result the residuals in the bottom panel do 
not exactly correspond to the differences between the data and models in the 
top panel.) Inset: likelihood functions for the mean of the residuals are shown 
for different subsets of the data, including the full sample, low- and high-mass 
subsamples, and the SFR divided subsamples. With one exception, the likelihood 
functions shown here include the galaxy 23.5+2-0.33, which drives some of the 
trends toward larger positive residuals (see text). The thin green curve shows 
the effect of removing this outlier. 

(A color version of this figure is available in the online journal.) 

the metallicities for our sample using the Maiolino et al. (2008) 
calibration of R23. 

Figure 7 compares our galaxies to the z — 0. 1 FMR reported 
in Mannucci et al. (201 1). On the top panel, we show the MZ 
relation, with our data points color-coded by three different SFR 
bins. The FMR, which is evaluated and shown for the same three 
SFR bins, is given by the relation (Equation (2) from Mannucci 
et al. 2011) 

12 + log(0/H) = 8.90 + 0.37m - 0.14s - 0.19m 2 

+0.l2ms — 0.054s 2 for p 032 ^ 9.5 (1) 

= 8.93 + 0.51(^0.32 — 10) for p 0.32 < 9.5, 

where m = log(M) — 10, s = log(SFR), and p 32 = log (M) — 
0.32 log(SFR). This exercise shows clearly that a planar relation 
exists to lower masses at z ~ 0.6-0. 7, with low SFR galaxies 
showing higher metallicities. To highlight the differences from 
the FMR, the bottom panel shows the residuals from the local 
plane. Here, we have included the uncertainty in the SED- 
derived SFR, as it enters into the metallicity that is predicted 


Sample 

Mean of the Residuals, fi 

Full sample 

0.016 + 0.020 

M > 10 9 ' 5 Mq 

-0.074 ± 0.083 

M < 10 9 ' 5 Mq 

0.046 ± 0.022 

-1.41 ^ log SFR < -0.42 

0.013 ± 0.035 

-0.36 ^ log SFR ^ 0.02 

0.086 ± 0.033 

0.08 ^ log SFR < 1.23 

-0.042 ± 0.029 

Without 23.5+2-0.33 

Full sample 

-0.009 + 0.021 

-0.36 ^ log SFR SC 0.02 

0.027 ± 0.043 

M < 10 9 ' 5 Mq 

0.018 + 0.025 


by the local FMR. Overall, we see that most galaxies have 
metallicities that fall within 0.2 dex of the local FMR. However, 
at low masses and low SFRs there are more objects with Alog 
(O/H) > 0, so it is reasonable to question whether these galaxies 
deviate slightly from the FMR. Therefore, we calculate the 
mean of the residuals, p = (Alog (O/H)), and its uncertainty 
by using a maximum likelihood estimation to account for 
the individual measurement errors. Under this approach, the 
probability of observing a galaxy with a residual v, and error a,- 
is pi = g-Cfi-AO / 2ct , Then, for the full sample the likelihood 
function becomes C(p) = J~[ p,. We evaluate this likelihood 
function between —0.4 < p < 0.4, taking the maximum and 
the 68% confidence intervals of dp) to represent the mean 
and its error. This estimate is made for the full sample, as well 
as low- and high-mass subsamples (divided at 10 9 5 M 0 ) and the 
subsamples in three bins of SFR. 

In the inset panel of Figure 7, we show the likelihood functions 
for the mean of the FMR residuals. These calculations indicate 
only insignificant deviations from the local FMR. In Table 4, 
we give the mean of the residuals and its uncertainty for the 
full sample and various subsamples. At first glance, we see two 
weakly significant deviations from the local plane: one among 
the lower-mass ( M < 10 9 5 M 0 ) galaxies, and the other for 
those with intermediate SFRs (green points). However, upon 
closer inspection we find that these trends arise solely because 
of one galaxy (23.5+2-0.33) with a metallicity 3.4a above 
the fundamental plane. (Furthermore, this galaxy has a strong 
influence on our calculation not because its residual is large, 
but rather because its errors are small.) Without this galaxy, the 
deviation disappears (see the thin green curve in the inset of 
Figure 7), and the mean of the residuals falls within 1.5a of 
zero for all of the subsamples listed in Table 4. Since any claims 
of a deviation from the FMR should be based on more than one 
object, we conclude that our data show no compelling evidence 
for evolution at 10 s < M/M 0 < 10 9,5 . This result lends greater 
leverage to the finding that the FMR is not evolving for higher- 
mass galaxies at these redshifts (Cresci et al. 2012). 

It is also interesting to examine the scatter about the FMR, 
as intrinsic scatter could indicate that there are additional 
physical processes that influence gas-phase metallicities. For 
intermediate-redshift galaxies, Cresci et al. (2012) report a 
scatter of 0.1 1 to 0.14 dex about the FMR. Yet, in our sample, 
the scatter appears to be larger at lower masses. We find an rms 
of 0.2 dex below M < 10 9 ' 5 Mq. Assigning some galaxies to the 
lower branch of R23 does not reduce this scatter. The amplitude 
of the residuals is larger (or about the same for a few galaxies) 
when the lower-metallicity solutions are adopted. In order to 
disentangle intrinsic scatter from increased measurement errors 
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in faint galaxies, we use a Monte Carlo simulation. First, we 
assume that there is no intrinsic scatter relative to the FMR. 
Then, focusing only on the galaxies with M < 10 9 ' 5 M q , 
we generate 10,000 mock realizations of the FMR residuals, 
perturbing the data in Figure 7 by their uncertainties. We 
calculate the rms about the FMR for each realization and find 
a mean of 0.23 dex, in good agreement with observations. This 
exercise implies that the amount of scatter that we observe can 
be entirely explained by our measurement errors. Therefore, we 
conclude that the intrinsic scatter about the FMR must be small. 
Repeating the simulation with additional intrinsic scatter allows 
us to place an upper limit on this scatter. We find that an intrinsic 
scatter of 0.16 dex produces an observed rms ^ 0.2 dex in only 
5% of Monte Carlo realizations. Remarkably, the scatter about 
the FMR for the present sample is significantly smaller than 
the 0.4 dex rms found for similar-mass galaxies in the local 
universe (Mannucci et al. 2011). We therefore conclude that the 
increased scatter in the low-mass, z = 0.1 FMR is likely the 
result of increased measurement uncertainties for faint galaxies. 

7. CONSTRAINING GAS FLOWS WITH THE 
MASS-METALLICITY AND SFR-M, RELATIONS 

As we introduced in Section 1 , the MZ relation is an important 
measure of galaxy evolution. Galaxies are not closed boxes; gas 
must be accreted from the IGM to fuel star formation (e.g., Erb 
2008; Bouche et al. 2010), and galactic outflows are commonly 
observed at both low and high redshift (e.g., Heckman et al. 
1990; Martin 2005; Henry et al. 2007; Soto et al. 2012; Kornei 
et al. 2012; Erb et al. 2012). The MZ relation provides an 
independent, albeit indirect, constraint on these gaseous flows. 
On one hand, the accretion of primordial gas will reduce the gas- 
phase metallicity in a galaxy. On the other hand, galactic winds 
can create a deviation from the closed-box model by lowering 
the gas fraction without changing the metallicity (Edmunds 
1990; Dalcanton 2007). Combining our data with the higher- 
mass MZ relation provides an important increase in dynamic 
range that allows for more stringent constraints on models. 

Our emission-line selection of galaxies does not strongly bias 
our sample or hamper our ability to test theoretical models. In 
Figure 5 we show that our sample is consistent with the DEEP2 
MZ relation reported by Zahid et al. (20 1 1 ; at least in the limited 
mass range where our data overlap). This agreement can be 
understood, since the FMR dependence on the SFR is relatively 
weak. For example, if the SFRs of our galaxies were biased 
high by a factor of several, then (following the Mannucci et al. 
2011 FMR) an “unbiased” sample would have metallicities that 
are approximately 0.1 dex higher. An increase in metallicity 
this large is implausible, as it would imply a lack of metallicity 
evolution inconsistent with other studies. Likewise, it is worth 
pointing out that the high signal-to-noise ratio (S/N) of the 
spectra shown in Figures 1 and 2 implies that our sample is not 
biased toward higher metallicities by a requirement for strong 
H/3 emission. The objects that were excluded due to low S/N all 
fell among the Class 3 objects discussed in Section 2.1; these are 
either extremely faint, low-mass galaxies with uncertain optical 
counterparts, or brighter galaxies that fell partially outside our 
slit. These galaxies have poorer quality spectra overall; it is 
unlikely that their exclusion biases our results. 

In order to further verify that our sample is not biased by 
emission-line selection, we place our galaxies on the SFR-M* 
relation (i.e., the star-forming main sequence) in the top panel 
of Figure 8. This comparison demonstrates that our galaxies 
do not have extreme SFRs for their stellar masses. Rather, 


log M ft /M s 



Figure 8. The SFR-M* correlation (top) and MZ relation (bottom) constrain 
models from Dave et al. (2012). Large red squares show the average mass, 
metallicity, and SFR for four mass bins: M > 10 9 ' 5 M 0 , 10 8 - 9 < M/Mq < 
10 9 ' 2 , 10 8 ' 5 < M/Mq < 10 8 ' 9 , and M < 10 8 ' 5 M e . Error bars on the binned 
data represent the standard deviation of the mean. We combine our data with 
binned higher-mass data from DEEP2: the z = 0.5-0. 7 SFR-M* correlation 
reported by Noeske et al. (2007; converted to a Chabrier 2003 IMF), and the 
Z ~ 0.8 MZ relation from Zahid et al. (201 1). MZ models are normalized by 
taking a nucleosynthetic yield that minimizes / 2 for the union of our data and 
the Zahid et al. data. SFR-M* models are not renormalized. The models are 
described in detail in Section 7.1. The halo mass scale shown at the top assumes 
the M*-Mhai 0 relation from Moster et al. (2010). 

(A color version of this figure is available in the online journal.) 

they lie on an extrapolation of the relation reported by Noeske 
et al. (2007). Our sample could only be significantly biased if 
the “true” SFR-M* relation turns over and becomes steeper at 
low stellar masses. While this trend is predicted to occur when 
photoionization heating suppresses star formation in low-mass 
halos, it is not expected to be important above Mhaio Z/ 10 9 M G 
(Okamoto et al. 2008). While beyond the scope of this paper, 
it is worth noting that we do not detect the flattening of the 
SFR-M* relation at low masses that is reported by Pirzkal et al. 
(2012). 

7.1. Theoretical Framework 

One set of models that have recently gained traction are those 
that adopt an equilibrium between mass inflow (Mj n ), outflow 
(M out ), and gas consumption via star formation (M*; Bouche 
et al. 2010; Dutton et al. 2010; Finlator & Dave 2008; Dave et al. 
2012). These models provide a relatively simple interpretation 
of the FMR as a manifestation of stochastic star formation. 
When a galaxy accretes gas, its metallicity is diluted and its 
star formation increases. The increased star formation serves to 
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consume the extra gas and return the galaxy to its equilibrium. 
Likewise, a pause in gas accretion will lower SFRs and increase 
gas-phase metallicities until accretion resumes. These models 
predict not only an MZ relation but also the SFR-M* relation 
that we have shown in Figure 8 above. Therefore, in the 
remainder of this section we leverage the two relations to gain 
the best possible constraints on the models. 

Quantitatively, we can express the equilibrium condition as 
(reproducing Equation (1) from Dave et al. 2012) 

— M out + AT. . (2) 

Following this assumption, Finlator & Dave (2008) and Dave 
et al. (2012) show that SFRs and metallicities can be written as 


4. y is the nucleosynthetic yield of oxygen, which is between 

0.008 < y < 0.021 (Finlator & Dave 2008). Since 
metallicity calibrations are uncertain (Kewley & Ellison 
2008), we take y as a free parameter as we search for a 
normalization that best matches our observed MZ relation. 
Hence, our comparison to the MZ relation concerns its 
slope, but not the normalization. 

5. 77 is the mass loading factor, which is defined as the ratio 
of the outflow rate to the SFR: q = M out /M*. In the 
comparison that follows, we adopt mass-loading factors 
that correspond (respectively) to momentum- and energy- 
driven winds (i.e., Murray et al. 2005; Finlator & Dave 
2008; Dave et al. 2011a, 2012; Peeples & Shankar 201 1): 


SFR = 


(" Mg rav 

(1 + rj)( 1 - a z ) 


(3) 


and 


_y L_ 

1 + rj 1 - a z 


(4) 


The quantities in Equations (3) and (4) are defined as follows. 


1. Zo is taken to represent the mass fraction of oxygen in 
the ISM. This quantity is converted to the units of 12+ 
log(0/H) by taking log(Z 0 ) = log(0/H) + log(3/4 x 
Mq/ M u ). Mo and M H are the atomic masses of oxygen 
and hydrogen. 

2. M grav is the cosmological baryonic accretion rate, taken 
from Dekel et al. (2009). Because this quantity (as well as 
others below) is expressed in terms of halo mass, we adopt 
the M*-Mhaio relation from Moster et al. (2010). 

3. £ is a quantity that represents preventive (rather than ejec- 
tive) feedback from gas heating. There are multiple forms 
of preventive feedback, and they combine multiplicatively. 
In the range of halo masses where our objects lie, the most 
important of these is gravitational heating from virial shock 
formation in accreted gas. Following Dave et al. (2012), we 
take the analytical form given by Faucher-Giguere et al. 
( 2011 ): 


Cgrav 



Hhalo 
10 12 Mq 


- 0.25 


(5) 


At high masses star formation quenching associated with 
supermassive black holes becomes important. While we do 
not aim to constrain quenching, we implement it following 
Dave et al. (2012): 


Cqucnch “ ( 1 + 1/3 (io' 2 U/ 0 ) ) ■ (6) 

The mass scale that we have adopted for £ quenc h is for z — 0, 
and while it may be higher at earlier cosmic times (Dekel 
et al. 2009), the exact scaling is unimportant in the present 
analysis. At the low-mass end, we take the photoionization 
feedback to be unimportant; as mentioned above, it is not 
expected to play a significant role at Mhaio > 10 9 M Q . 
Finally, Dave et al. (2012) do note that heating by galactic 
winds may be important, but the physics of this effect 
is poorly understood. As such, they adopt an arbitrary 
parameterization of e wm ds to consider its qualitative effects. 
As we will show below, this additional heating is not needed 
to reproduce the present data. 


i? = (M halo / 10 12 M o r * 1 / 3 

(7) 

>7 = (M halo / 10 12 M o r 2 / 3 . 

(8) 


The normalizations of these ij are taken from Dave et al. 
(2012). In addition to these two parameterizations of q, we 
compare to the case of no winds (q = 0), as well as q = 1, 
where the mass-loss rate is fixed to the SFR. 

6. Finally, the quantity az is the ratio of the metallicities 
of infalling gas and the ISM: Z m /Z ^ M . The inclusion of 
az allows for enriched inflows (possibly from previously 
ejected material now being re-accreted from the galaxy 
halo). Additionally, it appears in the model for SFR as 
a simple way to quantify the recycling of enriched halo 
gas into the ISM (i.e., an additional source of fuel for 
star formation). The parameterization of a z is taken from 
momentum-driven wind simulations in Dave et al. (201 la, 
2012 ): 


a z = (0.5-0.1z)(M*/10 10 M G )°' 25 (9) 

(although they note that it is a “crude” parameteriza- 
tion). While we aim to test models beyond those that are 
momentum-driven, a z is not reported for these simulations. 
Therefore we use this parameterization for all of the models 
with winds. For the no-wind case we consider a z = 0, since 
there can be no recycling of previously ejected material in 
this scenario. We consider the effects of modifying a z in 
the discussion below. 

7.2. Comparison of Equilibrium Models to 
Intermediate-redsliift Data 

Figure 8 shows how the equilibrium models compare to the 
intermediate-redshift MZ and SFR-M* relations. In order to get 
the best constraints on the models, we combine our sample 
with binned DEEP2 data representing the MZ relation for 
z ~ 0.8 galaxies (Zahid et al. 2011) and the SFR-M* relation 
for z ~ 0.5-0. 7 galaxies (Noeske et al. 2007). As we mentioned 
above, since metallicity calibrations and nucleosynthetic yields 
are uncertain, we fit for the yield that reproduces the MZ 
data. The yields that minimize x 2 (including uncertainties in 
both metallicity and mass 8 ) are y = 0.012, 0.013, 0.010, and 
0.008 for the momentum-driven, energy-driven, constant-wind 


8 We take X f = (OH, - /(M,)) 2 /^ , + /'(M ; ) 2 «t 2 _.), where f(M,) and 
f'(Mj) are the model and its derivative evaluated at M,-. yf is the sum over the 

i observations of yf . 
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(rj = 1), and no-wind (7/ = 0, a z — 0) models, respectively. We 
next describe how each model compares to the data. 

No wind. Not surprisingly, equilibrium models that exclude 
outflows are unable to reproduce the data. Most drastically, 
when rj — 0 and a z = 0, Equation (4) has no mass dependence. 
Hence, the predicted MZ relation is implausibly flat, and the 
slope of the SFR-M* relation is also too flat. Figure 8 shows 
that, without winds, low-mass galaxies make too many stars 
(and metals) and are unable to get rid of the metals that they 
make. It is worth noting that oe z is not strictly zero in the 
no-wind simulation (Dave et al. 2011b), since infalling gas can 
be enriched by other galaxies or could represent the accretion 
of a satellite galaxy. A non-zero a z would further increase the 
SFR and amplify the disagreement between the SFR-M* data 
and the no- wind model. 

Constant wind. A model where the mass-loss rate is equivalent 
to the SFR shows good agreement with the SFR-M* relation at 
all stellar masses. However, this model predicts an MZ relation 
that is too flat at low stellar masses. 

Momentum-driven wind. Of the models shown in Figure 8, the 
momentum-driven wind shows the best match to the MZ data. 
However, the metallicities predicted by this model turn up at 
higher masses, owing to the dependence on a z . Equation (4) 
shows that the predicted metallicities diverge as a z approaches 
unity at higher masses. While the disagreement between the 
model and data is only apparent in the highest-mass bin, the 
upturn remains implausible. Locally, the MZ relation flattens to 
higher masses (Tremonti et al. 2004), so the model shown in 
Figure 8 crosses the z ~ 0.1 relation at M ~ 10 10,5 M 0 . This 
observation suggests that the mass scaling for az given by Dave 
et al. (2012) may be too steep at high masses. 

In addition to the MZ relation, we compare the momentum- 
driven wind model to the SFR-M* relation in the top panel 
of Figure 8. In this case, the model slightly under-predicts the 
SFRs for the galaxies in the present study as well as for the 
lower-mass DEEP2 measurements. This discrepancy between 
the data and models at low masses is qualitatively similar to 
what is found at higher masses and lower redshifts. Dave et al. 
(201 lb) and Weinmann et al. (2012) find that both simulations 
and semi-analytic models are unable to reproduce the population 
of low-mass galaxies with high specific SFRs. 

Given the uncertainty in a z , we next consider whether 
modifying this variable can produce a momentum-driven wind 
model that is a good fit to both the MZ and SFR-M* data. 
If we remove the mass dependence of a z in Equation (9) 
(a, oc M°°), the predicted SFR-M* relation appears very 
similar to the constant-wind model (a good match to the data), 
but the predicted MZ relation (not shown) is far too flat. 
In summary, modifying a z does not improve the agreement 
between the data and the momentum-driven wind model. 
Energy-driven wind. The energy-driven wind models shown 
in Figure 8 are too steep compared to both the MZ data and 
the SFR-M* data. However, as we already noted, some of this 
steepening of the MZ relation occurs because a z may be too 
large at high masses. In fact, an energy-driven wind model 
with a z oc M® ° and y — 0.022 appears very similar to the 
momentum-driven wind model shown in both the top and bottom 
panels. For this model, the MZ relation model reduces to a 
straight line that runs through all of the data, and SFRs are 
slightly under-predicted at low masses. Ultimately, we cannot 
distinguish the momentum-driven model from an energy-driven 
model with a z oc M° °. However, we do note that this variable 


shows some mass dependence in the momentum-driven wind 
simulations (Dave et al. 2011a, 2011b), so it is unlikely that 
this dependence is non-existent for energy-driven winds. We 
therefore conclude that energy-driven winds are unlikely to 
explain the data. 

Recent simulations from Hopkins et al. (20 1 2) have suggested 
that outflows depend strongly on gas surface density, so that 
winds from low-density dwarf galaxies are dominated by energy 
driving from supernovae and stellar winds. More massive star- 
forming galaxies with higher gas densities, on the other hand, are 
predicted to prefer momentum-driven winds. The present data 
do not indicate a mass-dependent //; Figure 8 shows no evidence 
for steeper MZ and SFR-M* slopes toward lower masses. 

This analysis (and Figure 8) shows that there may be a tension 
between the data and the equilibrium model. On one hand, mod- 
els that best match the MZ relation under-predict the amount 
of star formation in low-mass galaxies (and even intermedi- 
ate masses probed by DEEP2). Likewise, models that best re- 
produce the SFR-M* relation over-predict the metallicities of 
galaxies. This trend can be understood as a generic property of 
the models, because both Equations (3) and (4) are proportional 
to ( 1 + rj)~ 1 (1 — a-) 1 ; changes in these model components alter 
the MZ and SFR-M* relation in the same manner. However, if 
preventive feedback from gas heating has been overestimated 
(C, underestimated in Equations (3) and (5)), then SFRs can 
be increased without modifying the MZ relation. In fact, Dave 
et al. (2012) note that £ glav is taken from simulations that do not 
include metal-line cooling (Faucher-Giguere et al. 2011), and 
raise the question of whether this effect is important. Alterna- 
tively, if C is correct, then the data may suggest a deviation from 
the equilibrium model. 

Ultimately, we conclude that this disagreement remains 
subtle. While we have argued that emission-line selection does 
not largely bias our results, we cannot rule out a small effect. 
As an example, if the SFRs presented in Figure 8 are high 
by 0.1 dex, we would expect an unbiased sample to have 
an SFR-M* relation that agrees better with the momentum- 
driven wind model at low masses. Additionally, according to 
the FMR, SFRs that are 0. 1 dex lower should be accompanied 
by metallicities that are around 0.03 dex higher, in even better 
agreement with the momentum-driven wind MZ relation. Even 
without emission-line selection bias, it is important to recall 
that the [On], H/J, and SED-based SFR indicators discussed in 
the Appendix differ systematically by up to 0.2 dex. However, 
modifying the SFR under these scenarios cannot satisfactorily 
bring the data in line with the momentum-driven model. Even 
in the higher-mass data alone, this model predicts a steeper 
SFR-M* relation than the DEEP2 result (Noeske et al. 2007); 
therefore, the discrepancy that we have identified is apparent 
in a magnitude-limited sample. We conclude that Ha follow- 
up spectroscopy of our sample (combined with a mass-selected 
control sample) can definitively assess the effects of emission- 
line selection, while simultaneously reducing errors on dust, 
SFR, and metallicity. We leave this work as the subject of a 
future paper. 

Finally, it is worth pointing out that systematic uncertainties 
in our metallicity calibration cannot explain the discrepancy be- 
tween the data and models. Kewley & Ellison (2008) show that 
local (SDSS) MZ relations have similar low-mass slopes under 
most calibrations. The Tremonti et al. (2004) calibration and 
the direct-method electron temperature calibration (Andrews & 
Martini 2013) are two exceptions — both produce MZ relations 
that are steeper than what we have used here. Hence, if we 
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calculated metallicities in the same way as Tremonti et al. or by 
using electron temperatures, the tension between the data and 
models would be amplified. On the other hand, the calibrations 
reported by Pilyugin (2001) and Pilyugin & Thuan (2005) give 
MZ relations that are flatter than most. However, we note that 
they also have metallicities around 0. 5-0.7 dex lower than those 
from KK04. These lower metallicities would introduce other 
difficulties, as they would require normalizing yields that fall 
below the plausible range indicated in Finlator & Dave (2008). 

8. CONCLUSIONS 

In this work we have placed the first constraints on the low- 
mass, intermediate-redshift MZ relation. By using emission- 
line selection in the COSMOS field, we are able to efficiently 
identify and measure the metallicities from 26 galaxies reaching 
masses of 10 s M Q at z ~ 0.65. Combined with previous 
measurements from magnitude-limited samples, these data 
extend our knowledge of the MZ relation to masses an order 
of magnitude smaller. This limit is comparable to the low- 
mass limit of the local MZ relation determined from SDSS 
data. Therefore, for the first time we are able to measure the 
metallicity evolution of intermediate-redshift galaxies at M < 
I 0 9 M q . Compared to the z ~ 0. 1 relation reported by Tremonti 
et al. (2004), we find an average metallicity that is 0.12 dex 
lower for galaxies with 10 85 < M/Mq < 10 90 . We show 
that this measurement is consistent with an MZ relation that 
evolves more strongly at low stellar masses. We interpret this 
mass-dependent evolution as consistent with downsizing trends, 
where higher-mass objects have less leverage to alter their gas- 
phase metallicities after most of their stars have been assembled. 

An important development in our understanding of metal- 
licity evolution is the discovery that the scatter in the relation 
can be reduced by accounting for star formation. This measure- 
ment was quantified as a plane, and dubbed the Fundamental 
Metallicity Relation (FMR) by Mannucci et al. (2010). Using 
the present sample, we find that a planar relation does indeed 
exist among low-mass, intermediate-redshift galaxies. Consis- 
tent with the findings of Mannucci et al., we see no evidence 
for evolution of the FMR, as our emission-line-selected sample 
falls in good agreement with their z ~ 0.1 FMR (Mannucci 
et al. 2011). Furthermore, while the lowest-mass galaxies in 
our sample show significant scatter with respect to the FMR, 
we attribute this scatter to measurement uncertainties for faint 
galaxies. Using a Monte Carlo simulation, we determine that 
the intrinsic scatter is consistent with zero and should be less 
than 0.16 dex (95% confidence) at M < 10 95 M 0 . 

The MZ relation is an important probe of galaxy evolution 
models, as accretion and galactic outflows modulate the gas- 
phase metallicities of galaxies. Hence, comparison to model 
predictions can help us to gain insight into galaxy formation. We 
have combined our MZ relation with the SFR-M* correlation 
that is also measured from our data, taking higher-mass data 
from the literature (Noeske et al. 2007; Zahid et al. 2011). 
We compare to the family of models outlined in Dave et al. 
(2012), where it is assumed that galaxies prefer to maintain 
an equilibrium between inflows, outflows, and star formation 
(Finlator & Dave 2008; Bouche et al. 2010; Dutton et al. 2010; 
Dave et al. 2012). We find that models that predict the MZ 
relation may under-predict the SFRs of low-mass galaxies, 
and at the same time, models that predict the SFRs of low- 
mass galaxies tend to over-predict their metallicities. While this 
finding could represent a breakdown of the equilibrium model 
in low-mass galaxies, it could alternatively be an indication 


that feedback from gas heating has been overestimated in 
simulations (i.e., Faucher-Giguere et al. 2011). To solidify this 
result, we have begun Ha and [N ii] spectroscopy to follow up 
the present sample and simultaneously measure a mass-selected 
control sample. These data will clarify the effects of emission- 
line selection and greatly reduce the statistical uncertainties in 
our metallicity, dust, and SFR measurements. 

To conclude, these observations have provided a valuable 
look at the metallicity evolution of low-mass galaxies outside 
the local universe. We look forward to the physical insights and 
new constraints that can be gained from larger samples and a 
better characterization of systematic uncertainties. 
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APPENDIX 

ASSESSING SFR MEASUREMENTS 

In Sections 6 and 7 we draw conclusions from SFRs that we 
derived from SED fits; therefore, it is important that we use 
reliable SFRs. In order to assess systematic uncertainties, we 
compare this measurement to two additional SFR indicators: 
H/J and [On]. For the H/J method, we use the dust extinction 
from the SED fits, scaled up by a factor of 2.3 to account for 
the difference between stellar and nebular extinction (Calzetti 
et al. 2000). We also correct the H/i fluxes for the small 
amount of stellar absorption, as discussed above. Next, the H/i 
luminosities are converted to Ha, assuming a Balmer decrement 
of 2.8. We then use the Ha-SFR calibration given by Kennicutt 
(1998), scaled appropriately to a Chabrier (2003) IMF. To derive 
SFRs from [Oil] luminosity, we use the calibration given by 
Moustakas et al. (2006), which includes an Mb dependence but 
does not require that the line measurement be corrected for dust. 

In order to make a direct comparison to the SED-derived 
SFRs, we must correct the emission-line fluxes for both slit 
losses and extraction aperture losses. Rather than model these 
losses (which would not account for uncertainties in the absolute 
flux calibration), we compare our data to narrowband imaging 
in COSMOS. First, we infer the emission-line flux from the 
narrowband (NB816) photometry, using the methods outlined 
in Ly et al. (2007) and Takahashi et al. (2007). Since the 
NB816 bandpass has a non-uniform throughput, we calculate 
a correction to this flux based on the relative throughput at the 
observed wavelengths of the emission lines. (This correction 
can be important, since it is not unusual that both of the [O m] 
doublet lines are found away from the central wavelength of 
the NB816 filter.) Because the COSMOS catalog contains total 
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Figure 9. Star formation rates derived from three different methods are compared. Blue points assume a constant star formation history, and the red points assume 
exponentially declining SFRs (r models) discussed in Section 3.2. The H/2-derived SFRs are corrected for dust extinction using these SED fits. (The center and right 
panels — where the SED-independent [O n] SFRs are compared — show that the r models tend toward slightly less dust.) The [O n] and H/3 SFRs are corrected by an 
average factor of 1.4 to account for slit losses. Systematic shifts between the different indicators are at most 0.2 dex. 

(A color version of this figure is available in the online journal.) 


magnitudes (under the reasonable assumption that the galaxies 
are unresolved with 1"5 FWHM spatial resolution), the line 
fluxes inferred from narrowband imaging are not subject to 
aperture losses. For 17 galaxies in our sample with emission 
lines that are bright enough to be detected in the narrowband 
imaging, we can infer the aperture losses by comparing these 
narrowband imaging fluxes to our spectroscopic fluxes. We find 
typical correction factors of 1. 1-2.0, with a median of 1.4. We 
adopt this average correction to explore systematic differences 
between star formation indicators. 

Figure 9 compares the SFRs derived from the SEDs, H/3, 
and [On]. In addition to the SED fits described in Section 3.2 
(which use exponentially declining star formation histories 
with r as a free parameter; red points), we also make the 
comparison for constant star-forming models (blue points). 
This exercise shows that the SFRs derived using different 
methods are correlated, so that we can easily distinguish the 
galaxies with high and low SFRs. However, Figure 9 also 
shows that there are systematic shifts as large as 0.2 dex 
between the different estimators. Determining the most accurate 
star formation indicator is beyond the scope of this paper. 
Ultimately, a systematic uncertainty of 0.2 dex does not affect 
our conclusions. In Section 6, we show that the scatter in 
the MZ relation is correlated with SFRs. As with the local 
FMR, metallicity depends only weakly on the SFR, so that 
a 0.2 dex offset in SFR translates to a 0.03 dex shift in the 
metallicity predicted by the local FMR (Mannucci et al. 2010, 
2011). This difference in metallicity is much smaller than the 
uncertainties on the FMR residuals shown in Figure 7. Likewise, 
in Section 7 we use our SFR- AT, relation to constrain models. 
However, 0.2 dex of systematic uncertainty is small compared 
to the dynamic range of SFRs shown in Figure 8. A systematic 
offset is unlikely to alter our conclusions, as the SFR- AT, trends 
discussed in Section 7 can be seen (albeit less obviously) from 
the higher-mass data alone. 
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